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The formation of Ps in our Galaxy is well measured, and has led to important and unanswered 
questions on the origin of the positrons. In principle it should be possible to form analogous 
systems from /r and r leptons, viz. true muonium and true tauonium. However the probability of 
formation for these systems is greatly reduced due to the intrinsically short lifetimes of the p and 
t leptons. Likewise, the decay of the atoms is hastened by the high probability of the constituent 
particles decaying. Nevertheless, if sufficient numbers of /r and r pairs are produced in high energy 
astrophysical environments there may be significant production of true muonium and true tauonium, 
despite the small probabilities. This paper addresses this possibility. 

We have calculated the pair production spectra of and r leptons from photon-photon anni¬ 
hilation and electron-positron annihilation in astrophysical environments. We have computed the 
cross-sections for radiative recombination and direct annihilation of the pairs, and the decay con¬ 
stants for the various allowable decays, and the wavelengths and energies of the recombination and 
annihilation signatures. In this way we have calculated the probabilities for the formation of true 
muonium and true tauonium, and the branching ratios for the various observable signatures. 

We have estimated the expected fluxes from accretion discs around microquasars and active 
galactic nuclei, and from interactions of jets with clouds and stars. We find that accretion discs 
around stellar mass black holes in our own Galaxy should have observable signatures at X-ray and 
7 -ray energies that are in principle observable with current observatories. 

PACS numbers: 34.80.Lx, 36.10.Dr, 36.10-k, 95.30.Cq, 95.30.Ky, 12.20.Ds, 13.66.De 


I. INTRODUCTION 

In the years 1928 1931 Dirac introduced the con¬ 

cept of antimatter [l]-|3|. In 1932 Anderson showed that 
these theoretical concepts were real with the detection 
of the positronQ. Then in 1934 Mohorovicic published 
a paper on the possibility of observing positronium (Ps) 
from astrophysical sources via its recombination lines[l . 
This prescient paper went largely unnoticed, and the ex¬ 
istence of Ps was predicted independently and separately 
by Pirenne RuarkjT] , Landau (unpublished work re¬ 
ferred to in Alichanian and Asatiani Q) and Wheeler [|E|. 

To date, the observational challenge laid down by Mo¬ 
horovicic, viz. to detect astrophysical sources of Ps from 
its recombination lines, has not been met. However the 
possibility of detection of Ps recombination lines has been 
raised on subsequent occasions inspired by advances in 
technology [10j-ll2j . Despite the non-detection of Ps re¬ 
combination lines, there is clear evidence of its existence 
in the Galaxy from its annihilation signature[l3j: see the 
recent review by Prantzos et al. |l4j. 

The observations of annihilating electrons and 
positrons in our Galaxy have raised several important 
questions for astrophysics, the foremost being what are 
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the sources of the positrons? Can they be explained with 
conventional astrophysics (e.g. high energy sources) or 
are more exotic sources required (e.g. annihilating dark 
matter)? Why is the annihilation radiation centered so 
strongly on the nucleus of the Galaxy with relatively lit¬ 
tle disc emission? How far do positrons disperse through 
the Galaxy before annihilation? 

In this paper we address the question of what other ob¬ 
servable signatures there are, arising from other leptonic 
atoms. For example, it is possible in principle to form 
atoms of fi~ - fj, + and r~ r + . However, in these cases 
the observable signatures are complicated by the fact 
that the constituent particles are themselves unstable, 
and may decay prior to annihilation or transitionsfl5|. 
It is important to note that even improbable signatures 
may be significant if the rate of production of the atom 
is high enough; this situation is often very relevant in 
astrophysical situations. For example the hydrogen 21 
cm emission line that results from the spin-flip transition 
between the hyperfine splitting of the 2 S ground state 
has been hugely exploited in radio astronomy since the 
1950sQSU7|, although it has a transition probability of 
only 2.9 x 10 -15 s -1 . 

Even faint signatures would be of interest. First, they 
would be the first such observations of /r - - or t~ 
t + , either in astrophysics or in the laboratory. Secondly 
they would place constraints on the energetics and rate 
of production of /r or r in astrophysical sources. Thirdly, 
the rate of - fi + or t~ t + production, compared 
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with estimates of /x or r production rates, would place 
constraints on the physical conditions during recombina¬ 
tion. 

Such measurements could be very useful. Taking the 
example of Ps, the rate, distribution and energetics of for¬ 
mation are all well measured, but the identity of positron 
sources still remains unclear. This is largely because it 
is difficult to reproduce the distribution of Ps emission, 
and also because there are many e + candidates, which 
are difficult to distinguish observationally. Observations, 
or even non-detection, of p~ - /x + or t~ ■ r + could place 
further constraints on the competing models for the ori¬ 
gin of Galactic positrons. 

The study of such processes is important, even if they 
are unlikely. Anomalous spectral signatures are routinely 
discovered, for example the recent detection of an uniden¬ 
tified emission line in the X-ray spectra of clusters of 
galaxies at 3.5 keVjjjl, [lsj]. We note that this line is 
not attributable to p~ - p + or t~ t + , although these 
species are expected to produce X-ray emission lines, as 
we show. Nevertheless, in such cases it is important to 
rule out prosaic, but poorly understood processes before 
accepting more exotic physics. The formation of p~ 
/x + and t~ - r + in astrophysical environments has so far 
received very little attention, but should not be rejected 
a priori on the basis of the short lifetimes of the /x and 
r. 

The production of p~ - p + and t~ t + in electron- 
positron colliders is difficult due to the rapid decay times 
of the particles, although methods have been proposed 
for the creation of p~ - /x + [2(|. It may be possible to by¬ 
pass these difficulties if large numbers of leptonium atoms 
are produced in high energy astrophysical environments. 
This is an important possibility since neither p~ - /x + nor 
t~ - t + have ever been observed, and together they are 
the most compact pure QED systems that could exist. 
An astrophysical detection would be more difficult to in¬ 
terpret than a collider detection, but would nevertheless 
constitute an important basic test of QED. 

We wish to make a first step in addressing this chal¬ 
lenge by calculating the expected observable signatures 
of astrophysical sources. Consider a source of particle - 
anti-particle lepton pairs. Their possible fates are, (i) the 
particle immediately decays, if it is or t±, (ii) it di¬ 
rectly annihilates with its antiparticle, (iii) it radiatively 
recombines with a free particle, then annihilates, (iv) it 
forms a new atom via charge-exchange, then annihilates. 

Ps clearly has the highest probability of forming, 
since its constituent particles are stable, and can travel 
through the interstellar medium until they form Ps 
through charge-exchange or radiative recombination, al¬ 
though they may also directly annihilate in-flight or when 
thermalized. The probabilities for Ps formation have 
been discussed in detail[ijJ ( 2 ]], |22| , and Ps formation 
at the Galactic centre is confirmed by the observation of 
the 511 keV annihilation signature and the three-photon 
triplet annihilation continuum[23]. from which a Ps for¬ 
mation fraction of ss 97 per cent is inferred. See Prantzos 


et al. M for a recent review of Ps astrophysics. 

The formation of any atom having a /i or r lepton is 
much less likely, since the decay times are so short that 
no travel of the particles is possible. Thus, we only con¬ 
sider the formation of true muonium and true tauonium 
(hereafter referred to as M and T01) in situ, i.e. the par¬ 
ticle - anti-particle pairs must be created with an energy 
less than the ionisation energy, so that the particle forms 
immediately. A summary of the basic properties of M, 
T, and Ps, to which frequent comparisons are made, is 
given in Tabled] the calculations or references for these 
values can be found throughout the paper, in the sections 
given in the final column of the table. 

To calculate the likelihood of forming M and T, we 
begin by reviewing the possible mechanisms for pair- 
production which are relevant in astrophysical environ¬ 
ments (section m ■ We then use the relevant cross- 
sections for pair-production, to calculate the pair- 
production spectra, and from these calculate the fraction 
of pairs which have sufficiently low energy to form an 
oniurn immediately, i.e. the fraction of pairs produced 
with kinetic energy less than the ionisation energy of 
the atom ('section IlIIl) . These pairs may either directly 
annihilate, or radiatively recombine, and we calculate 
the branching fraction for these processes in section HVl 
For those that do radiatively recombine, we calculate the 
cross-sections for radiative recombination, and hence the 
fractions which recombine into the nLth quantum level in 
section ITV Al Then in section[V]we calculate the different 
decay channels and in section IVll the resulting branching 
ratios for an onium in the nth level, and thus the total 
branching fractions for specific observational signatures. 
The signatures themselves are discussed in section IVIII 
This brings us to a position where, for a given rate of 
pair-production, we can estimate the luminosity of each 
of the different decay channels, and we apply these cal¬ 
culations to specific astrophysical sources in section lVIIIl 
We give our conclusions in section ||X] 


II. PAIR-PRODUCTION PROCESSES IN 
ASTROPHYSICAL ENVIRONMENTS 

The formation of leptonium other than Ps is not com¬ 
plicated by possibilities of escape, in-flight formation, 
charge-exchange etc.; the lifetimes of the particles them¬ 
selves are so short (§0, that if an atom is not formed 
immediately after the creation of a particle, then the par¬ 
ticle itself will decay before any of these processes can 
occur. Ultra-relativistic particles may live longer due to 
time dilation, but as for e + [24j , the annihilation signature 
will give rise to very weak continuum emission. 

We consider the following production processes for lep- 
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TABLE I: Summary of the main properties of Ps, M and T. 

Ps M 


T 


Section 


Rest mass/annihilation energy (MeV) 
Ionisation energy (eV) 

Bohr radius (m) 

Decay time of constituent particles (s) 
77 Annihilation time (ground state) (s) 

e“e + Annihilation time (ground state) 
Recombination time (nL —> n'L') (s) 


Recombination energies (keV) 




0.511 



105.66 

1406.6 




6.8 



1784.1 

23751.4 

m 



1.058 

x 10 -10 

5.199 

x 10 -13 

3.044 x 10~ 14 

Nvc| 



OO 



2.197 

x 10 -6 

2.874 x 10~ 13 

Nvb| 

Singlet 


1.2 x 

10 

r 10 

6.0 x 

HT 13 

3.6 x 10 -14 

MVA| 

Triplet 


1.4 x 

10 

r 7 

6.7 x 

10” 10 

4.0 x 10~ n 

MVA| 

Triplet 


- 



1.8 x 

10~ 12 

1.1 x 10~ 13 

MVA| 

2 1 —t 1 0 

(Lyman 

a) 3.191 

X 

10 -9 

1.543 

x 10 -11 

9.18 x 10 -13 

Mv c| 

31-110 

(Lyman 

/3) 1.195 

X 

10 -8 

5.780 

x 10 -11 

3.44 x 10~ 12 

Mv c| 

31^20 


8.905 

X 

10 -8 

4.307 

x 10 _1 ° 

2.56 x 10” 11 

Mv c| 

30->2 1 


3.166 

X 

10" 7 

1.531 

x 10 -9 

9.11 x 10” 11 

Mvc| 

3 2 -s- 2 1 


3.092 

X 

10 -8 

1.495 

x 10 _1 ° 

8.89 x 10 -12 

Mv c| 

4 1 -> 2 0 


2.068 

X 

10- 7 

9.999 

x 10 _1 ° 

5.95 x 10" 11 

Mv c| 

4 0 -> 2 1 


7.753 

X 

10- 7 

3.750 

x 10 -9 

2.23 x 10 _1 ° 

Mv c| 

4 2 -> 2 1 


9.692 

X 

10 -8 

4.687 

x 10 _1 ° 

2.79 x 10 -11 

Mv c| 

Lyman a 


5.102 

X 

10 -3 

1.055 


17.7 

Mv c| 

Lyman /3 


6.047 

X 

10 -3 

1.250 


21.0 

Mv c| 

B aimer a 


9.448 

X 

10 -4 

0.195 


3.3 

lvc| 

Balmer /3 


1.276 

X 

10 -3 

0.264 


4.4 

! JVC| 


tons, l, 


7 r — y l + z//, 

(i) 

7T - *” — y l^~ 

( 2 ) 

7 + 7 —>■ l~ + / + , 

(3) 

“+e + -> /"+/+. 

(4) 


The pion decay processes will occur in collisions of highly 
energetic (> 200 MeV) cosmic rays with protons in the 
ISM[25j. The third process requires very high photon en¬ 
ergies and will occur in the high energy environments sur¬ 
rounding compact objects, such as black-holes and neu¬ 
tron stars[25|. The last process may occur if the e + have 
large kinetic energy, and collide with e~ in the ISM, e.g. 
if an e 1 * 1 jet from an AGN or microquasar collides with 
an interstellar cloud of gas or a star. 

We can immediately dismiss the pion decay processes 
as irrelevant for r leptons, since the r are more mas¬ 
sive than the Tt^. We can also dismiss the pion decay 
processes as irrelevant for the formation of M, on two 
grounds. First, only a single particle is produced, and 
the probability of radiative recombination with another 
free n before it decays is negligible. Secondly, the en¬ 
ergetics are unfavorable. Since the reaction involves a 
single particle of known mass decaying into two parti¬ 
cles, the energetics are precisely determined in the zero- 
momentum frame; the p, is created with a kinetic energy 
of 4.12 MeV, which is much higher than the ionisation en¬ 
ergy of M, ss 1.4 keV. Thus the created /i must first lose 
energy, then collide with its antiparticle, before either 
particle decays, the probability of which is negligible. 

This leaves us with photon-photon pair production and 
electron-positron pair production. In both these pro¬ 
cesses the p, (or r) will be created with typical ener¬ 
gies much higher than the ionisation energy of M (or 
T) (» 1.4 keV for M, and ft 23.7 keV for T). Only the 


small fraction of pairs whose total kinetic energy in the 
zero momentum (z.m.) frame is less than the ionisation 
energy can form an onium, i.e., 


T1+T2 < £jon, 

M 4 


^32%* 


(5) 

( 6 ) 


where the right hand side is the ionisation energy in S.I. 
units, m is the mass of the lepton, /i is the reduced mass 
of the atom, q is the electron charge, and 7 is the Lorentz 
factor of the produced pairs in the z.m. frame. Therefore 
the maximum Lorentz factor in the z.m. frame that the 
particles can have and still form an onium is, 


7lim — 1 + 


1287r 2 c 2 e^ 2 


= 1 + 6.656 x 10 


-6 


(7) 


Thus, only very low energy pairs can produce an onium, 
and 7 n m is independent of mass. We now estimate the 
fraction of pairs produced with 7 < 7 ii m , for various for¬ 
mation mechanisms. 


III. THE PAIR-PRODUCTION SPECTRUM 

Throughout this section we will refer to /r and M, but 
the arguments are identical for r and T. 


A. Pair production by photon-photon annihilation 

For two photons with energy Ei and E 2 colliding at 
an angle 0 , pair production will occur if, 


E\E 2 7 


2 ( mc 2 y 
1 — cos t 


( 8 ) 
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where m = 2/x is the mass of an individual particle. The 
total energy of either produced particle, E 0 is given by, 

El = 7 2 m 2 c 4 = ^^(l-cosg) ^ (9) 

where 7 is the Lorentz factor of either of the created 
particles in the z.nr. frame. 

Many authors (26l - l30j have evaluated the pair- 
production spectrum (usually of e 41 ) based on various 
assumptions about the distributions of photon energies. 
Here we repeat these calculations, except we integrate the 
resulting spectra in order to obtain the fraction of pro¬ 
duced pairs which have sufficiently low energy to form 
M, i.e. those which obey equation [3 

The total cross-section for the reaction 7+7 —> 
is 0 ], 


TABLE II: The fraction of pairs produced via photon-photon 
annihilation with sufficiently low energy to form an onium, 
for photon distributions with a power-law energy density with 
index —a. 

a f 1 ~ l 

1 2.0 x 1CT Y 

1.5 3.8 x 10" 7 

2 6.1 x 10 -7 

2.5 8.9 x 10 -7 

“Since the fraction / 77 is independent of mass this is the fraction 
for any onium, Ps, M or T. 

to specify only the produced pairs which satisfy equa¬ 
tion [TJ or 

f^lmax = Fj loo = OO (15) 


er 


F _ - _ 

c 327reQm 2 c 4 

(1 - p 2 ) ((3 - /3 4 ) In - 2/3(2 - /3 2 )) (10) 


F r = 


1 - e 


( 11 ) 


where P = u/c, and u is the velocity of the outgo¬ 
ing particles in z.m. frame, and F c is the Sommerfeld- 
Sakharov correction to the cross-section, due to the at¬ 
tractive Coulomb force between the created /i , which 
applies near the threshold energy [32| - l34| . Now consider 
a source of photons whose number density per unit en¬ 
ergy is given by ni(Ei). Let these photons interact with 
an isotropic photon gas whose number density per unit 
energy is given by n 2 {E 2 ). The fraction of photons in the 
isotropic gas which lie within a differential cone at angle 
9 and width d9 to the incoming photon is 1 sin0dd. The 
reaction rate in the lab frame, i.e. the number of pairs 
created per unit volume per unit time is thus given by, 




r^in 


1 " ^lmin 


n\{Ei)n 2 (E 2 )<j(P)-(l - cos 9) 
AE\ d cos 9 dE 2 


( 12 ) 


where the factor c(l — cos 9) is the relative velocity of the 
incoming photons, and the lower limit on the integral 
over E \, 


E 


lmin — 


2 m 2 c 4 

E 2 {\ — cos 9) 


(13) 


from equation]®] ensures pair production is possible. The 
upper limit on E\ is either 


L'lmax — Eim 


2 7iL w2c4 

£ 2(1 — cos 9) 


(14) 


for all produced pairs. The fraction of pairs which can 
produce M, / 77 , is thus given by the ratio of R integrated 
over both these limits. 

We have evaluated equation [ 12 ] numerically for differ¬ 
ent cases of photon energy distribution. Because we are 
interested in the ratio of R for different limits in 7 , and 
because Ei and E 2 are related to 7 via equation EH the 
form for the second distribution n(E 2 ) does not matter; 
integrating over E 2 leads to a different constant which 
cancels when taking the ratio. Thus, / 77 is identical for 
a specific power-law interacting with any isotropic pho¬ 
ton gas. The fraction / 77 for different power-law indices 
are given in Table [TT] 


B. Pair production by electron-positron 
annihilation 

The total cross-section for the process e~ +e + —>• /i _ + 
/i + is given by [H|, 


na 2 H 2 p (P 4 - 4/3 2 + 3) 

° = F * - 6 MV- (16) 

where M is the muon mass, and F c is the Sommerfeld- 
Sakharov factor given by equation [lT] The fraction of 
those muons produced which have sufficiently low energy 
that they could immediately form M now given by 


R 


ee 



ni(Ei)n 2 (E 2 )a(p) 1 - 
dEi dcos0 dE 2 


(17) 


where v is the relative speed of the electron and positron 
in the lab frame and ni(Ei) and n 2 (E 2 ) are now the 
distributions of the lab frame number density per unit 
energy of the electrons and positrons. 

The lab frame energies Ei and E 2 are related to the 
Lorentz factors of the produced pairs by equation IA.31 
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TABLE III: The fraction of pairs produced via electron- 
positron annihilation with sufficiently low energy to form an 
onium, for electron distributions with a power-law energy den¬ 
sity with index —a. 

° _ fee 

1 5.7 x 10“' 

1.5 9.1 x 10 -7 

2 1.3 x 10 -6 

2.5 1.7 x 10 -6 

“Since the fraction f ee is independent of mass this is the fraction 
for any onium, Ps, M or T. 


and the relative velocity in the lab frame, v, is given by 
equation IA.4I Therefore changing variables to integrate 
over 7 instead of Ey, we can calculate i? e e 7r in the limit 
7 max = 7iim, and R eeoo in the limit 7 max = oo, and thus 
the fraction of produced pairs with sufficiently low energy 
to form M is given by, 


fee 


Re 


Re 


(18) 



FIG. 1: The fraction of thermalised pairs which will either ra- 
diatively recombine (R.R.) or directly annihilate (D.A.), ne¬ 
glecting charge-exchange, as a function of temperature, for e 1 * 1 
(continuous line), g ± (dashed line), and r ± (dotted line). 


As for photon-photon annihilation the exact energy den¬ 
sity of E 2 does not matter, since it cancels when taking 
the ratio, and the result is independent of the mass of 
the produced particles. We have calculated f ee for vari¬ 
ous power-law indices, where we evaluated the integrals 
numerically, and we show the results in Table IIIII 


IV. CROSS-SECTIONS FOR RADIATIVE 
RECOMBINATION AND DIRECT 
ANNIHILATION 

In the previous section we calculated the fraction of 
pairs produced, either via photon-photon annihilation, or 
by electron positron annihilation, which have sufficiently 
low energy to form Ps, M or T immediately. However, 
the fraction of these pairs which will form an onium is 
less still. As for PsyUJ there is a still a probability that 
some of these pairs will directly annihilate; though in the 
case of g and r we can neglect charge exchange since 
the probability of the particles decaying before meeting 
a neutral atom is large. 

We have calculated the total reaction rates for ra¬ 
diative recombination and direct annihilation following 
Gould [ 2 H (see also Gould ‘36J). Figure [Q shows the 
fraction of pairs which radiatively recombine, or di¬ 
rectly annihilate, g ^ a , as a function of temperature. For 
e^ radiative recombination dominates below ft 10 6 K, 
and for g and r, radiative recombination dominates at 
all astrophysically relevant temperatures. 


A. Cross-sections for radiative recombination onto 
the nLth level 

For those pairs which do radiatively recombine to form 
an onium, we require the fraction which recombine into 
the nLth level, such that we may then determine the sub¬ 
sequent decay. Wallyn et al. [22j give the cross-sections 
for the formation of Ps via radiative recombination as a 
function of the quantum levels n and L and the relative 
energy of of the electron and positron. We have repeated 
these calculations for Ps, M and T, and show the results 
in Figure [2] Thus we may calculate the fraction of oni- 
ums which form in the nLth energy level via radiative 
recombination as 


71=00 jL=71 — 1 

E E a ( nL ) 

n=l L —0 

In practice we cannot sum to n = 00 , and instead choose 
a suitably high value of n = nn m , such that the fraction 
fvnL asymptotes to a constant value. We computed f arlL 
for various limits, ?zn m , and found that ny lm = 300, 500 or 
1000 is sufficient for Ps, M, or T respectively. Using these 
limits we have calculated the fractions f (JnL for Ps, M and 
T as a function of temperature, and give the results for 
the first six energy levels in Table IIVI 

V. DECAY CHANNELS 

We have now calculated the fraction of pairs pro¬ 
duced via photon-photon (§ IIII AD . or electron-positron 
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TABLE IV: The fraction of those atoms which radiatively recombine to form Ps, M or T, which do so in the level nL. 


n 

L 

10 3 

10 4 

Ps 

10 s 

10 6 

10 7 

10 3 

Temperature (K) 
M 

10 4 10 5 10 6 

10 7 

10 3 

10 4 

T 

10 s 

10 6 

10 7 

1 

0 

0.296 

0.459 

0.702 

0.814 

0.830 

0.153 

0.193 

0.263 

0.396 

0.628 

0.123 

0.145 

0.182 

0.244 

0.358 

2 

0 

0.043 

0.068 

0.100 

0.105 

0.104 

0.022 

0.028 

0.039 

0.059 

0.093 

0.018 

0.021 

0.027 

0.036 

0.053 

2 

1 

0.116 

0.143 

0.068 

0.009 

0.001 

0.062 

0.078 

0.105 

0.140 

0.102 

0.050 

0.058 

0.073 

0.097 

0.133 

3 

0 

0.015 

0.023 

0.031 

0.031 

0.031 

0.008 

0.010 

0.013 

0.020 

0.029 

0.006 

0.007 

0.009 

0.012 

0.018 

3 

1 

0.044 

0.055 

0.025 

0.003 

0.000 

0.023 

0.030 

0.040 

0.053 

0.038 

0.019 

0.022 

0.028 

0.037 

0.051 

3 

2 

0.048 

0.037 

0.003 

0.000 

0.000 

0.027 

0.034 0.045 

0.046 

0.010 

0.022 

0.026 

0.032 

0.042 

0.049 

4 

0 

0.007 

0.010 

0.013 

0.013 

0.013 

0.004 

0.004 0.006 

0.009 

0.013 

0.003 

0.003 

0.004 

0.006 

0.008 

4 

1 

0.021 

0.026 

0.011 

0.001 

0.000 

0.011 

0.014 0.019 

0.025 

0.017 

0.009 

0.011 

0.013 

0.018 

0.024 

4 

2 

0.030 

0.023 

0.002 

0.000 

0.000 

0.017 

0.021 

0.028 

0.029 

0.006 

0.014 

0.016 

0.020 

0.026 

0.030 

4 

3 

0.020 

0.007 

0.000 

0.000 

0.000 

0.013 

0.016 

0.020 

0.013 

0.001 

0.010 

0.012 

0.015 

0.019 

0.017 

5 

0 

0.004 

0.006 

0.007 

0.007 

0.007 

0.002 

0.003 

0.003 

0.005 

0.007 

0.002 

0.002 

0.002 

0.003 

0.005 

5 

1 

0.012 

0.014 

0.006 

0.001 

0.000 

0.006 

0.008 

0.011 

0.014 

0.009 

0.005 

0.006 

0.007 

0.010 

0.013 

5 

2 

0.018 

0.014 

0.001 

0.000 

0.000 

0.010 

0.013 

0.017 

0.017 

0.004 

0.008 

0.010 

0.012 

0.016 

0.019 

5 

3 

0.017 

0.007 

0.000 

0.000 

0.000 

0.011 

0.014 0.017 

0.012 

0.001 

0.009 

0.011 

0.013 

0.017 

0.015 

5 

4 

0.008 

0.001 

0.000 

0.000 

0.000 

0.006 

0.008 

0.009 

0.003 

0.000 

0.005 

0.006 

0.007 

0.009 

0.005 

6 

0 

0.002 

0.003 

0.004 

0.004 

0.004 

0.001 

0.002 

0.002 

0.003 

0.004 

0.001 

0.001 

0.001 

0.002 

0.003 

6 

1 

0.007 

0.008 

0.003 

0.000 

0.000 

0.004 

0.005 

0.007 

0.008 

0.005 

0.003 

0.004 

0.005 

0.006 

0.008 

6 

2 

0.012 

0.009 

0.001 

0.000 

0.000 

0.007 

0.008 

0.011 

0.011 

0.002 

0.005 

0.006 

0.008 

0.010 

0.012 

6 

3 

0.013 

0.005 

0.000 

0.000 

0.000 

0.008 

0.010 

0.013 

0.009 

0.000 

0.007 

0.008 

0.010 

0.012 

0.011 

6 

4 

0.009 

0.001 

0.000 

0.000 

0.000 

0.007 

0.009 

0.010 

0.004 

0.000 

0.006 

0.007 

0.008 

0.010 

0.006 

6 

5 

0.003 

0.000 

0.000 

0.000 

0.000 

0.003 

0.004 0.004 

0.001 

0.000 

0.002 

0.003 

0.004 

0.004 

0.001 


(j ump . annihilation which have sufficiently low energy 
to form Ps, M or T, and of those which do, the frac¬ 
tion which radiatively recombine or directly annihilate 
(§hyd, and of those which radiatively recombine to form 
an oniurn, the fraction which form in the nLth level 

(§UYU) ; 

An onium in the nLth level may decay via annihilation 
of its constituent particles, if L = 0; it may radiatively 
transition to a lower quantum level if n > 1; or if it is M or 
T, either one of the individual constituent particles may 
decay from any level. We wish to calculate the branching 
fractions for these various decay channels, and we begin 
by calculating the lifetimes for each process. 


conservation of momentum. However, for M and T the 
triplet state can also decay via one photon, which then 
decays into e 1 * 1 pairs, i.e. M —> 7 * —>- e + e“. 

The decay rate of an onium in the nth level of the 
singlet state is given by (39], 


Ti = 


1 a 5 nc 2 


( 21 ) 


where /i is the reduced mass of the onium. 

For triplet states the lowest order decay rate is jdOlUl 


12 ,, 

r3 = —^ 
n 6 97r 


2/ic 2 


( 22 ) 


A. Annihilation 


For M and T the decay of the triplet state into electron- 
positron pairs has a decay rate given by®, 


Invariance under charge-conjugation leads to the fol¬ 
lowing selection rule for the decay of an onium into n 
photons, 


(-l) i+s = (-1)™, (20) 


a 5 11c 2 
SHn 3 


(23) 


B. Decay 


where l and s are the orbital angular momentum and 
spin quantum numbers respectively [13, [H]. Further¬ 
more, since the electron and positron (or /i~ and ii + 
etc.) only overlap in the L = 0 state, annihilation into 
photons is only possible (barring negligible higher order 
decay processes) from the singlet 1 Sq state or from the 
triplet 3 Si state. The 1 Sq will decay into an even num¬ 
ber of photons, with 2 being the most probable. The 
3 Si state will decay into an odd number of photons, with 
3 being the most probable, since 1 is forbidden due to 


Both the /i and the r lepton are intrinsically unstable 
and will decay to lighter particles. We take the decay 
rates to be, T ^ = 1/T ,,. = 2.197 x 10 -6 s and T t = 
l/r T = 2.874 x 10 -13 s [H]. 

C. Radiation 

The coefficient for a radiative transition of a hydro- 
genic atom from the state n'L' to a lower state nL, 
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Energy (eV) 

(a)The total cross section for levels nL 



Energy (eV) 

(b)The total cross section for level n, summed over all the 
L sub-levels. 


A n 'L',nL is given byjHHIH 


An'L',nL = 4ttV 


87raao\ 

Ik 5 "" ) 


Ma x(L, L ') 
(2 L> + 1) 


\p(n'L',nL )\ 2 , 
(24) 


where v is the frequency of the transition, 


, 1 1 
v = cRea I ,7-~ 


(25) 


77, 

Roo is the Rydberg constant for the particular onium, 


(Me 4 ) 

8£q ch 3 


(26) 


and ao is the Bohr radius of the onium, 


do 


4:Treoh 2 

9 ' 

fie z 


(27) 


FIG. 2: (a ) The cross section for the radiative recombina¬ 
tion of Ps (black), M (red) and T (blue), onto the quantum 
levels nL for n = 1 to 3 as a function of the relative energy 
of the combining particles. The lines are coded (1,0) thick, 
(2,0) normal, (2,1) thin, (3,0) dashed, (3,1) dotted, (3,2) dot- 
dashed. (b) The total cross sections for level n, from n = 1 
to 10, summed over all the L sub-levels. 


The dipole matrix elements are 


given 


by0, 


\p{n'L — l,ni| 2 = 


(-!) n/ 1 / (n + L)\{n' + L- 1)! (4 nn') L+1 _ , n+n '_ 2L -2 

4(2 L - 1)! V (n-L- l)!(n' - L)\ (n + n ') n + n ‘ 1 ] 


X 2^1 


— 7i T L T 1, — 7}! T T, 2T, 


—4 Tin' 


(n — n') 2 


71 — 71 

7i + n' 


2 -Pi 


—7i + L — 1, —7i' + L, 2 L, 


—4 nn' 


(■n — n ') 2 


(28) 























which is correct when L > L'\ if L < L' then nL and 
n'L' are swapped. These can be adapted for any oniurn 
simply by replacing the reduced mass, /r[22j. 

The radiative coefficient for a transition from a level 
n' L' to any lower level is thus given by, 


Similarly the probability for leptonium in the level n'L' 
to decay via two-photon annihilation is, 


= 

' L '^ 4 A n , v . 


(35) 


n=n' — 1 

A n 'L’ = ^ ^ A n : L :, rlL . 

n—1 L=L'± 1 


VI. BRANCHING RATIOS 


(29) 


Likewise for three photon annihilation, 


Pn 


3 S(L') r 3 


T / =- 

in A A , T / 


(36) 


and annihilation into an electron-positron pair 


A. Branching ratios for a particular level, nL 

Consider leptonium in the level nL. It may decay via 
annihilation (two photon, three photon or eA , depending 
on the state and the type of atom), radiative transition, 
or the decay of either constituent particle (for M and 
T). The decay rates for these channels were calculated 
in the previous section. Thus we may now calculate the 
probabilities for each process. 

Recalling that annihilation can only take place from 
the states with L = 0, and also that the annihilation 
mechanisms differ for singlet states, which decay into two 
photons, and triplet states, which decay either into three 
photons, or into e± pairs via a single photon, the total 
decay constant for leptonium in a level nL via any decay 
process is given by, 

^-niaing = A n L + 5(L)T 1 +2T^ t , 1 5' o (30) 

AnLtrip — A nL +S(L)T 3 + 5(L)r ee + 2T^ T , 3 $i(31) 

where 5(L) is a function which is 1 when L = 0, and zero 
otherwise. 

Therefore the probability of making a particular ra¬ 
diative transition from a level n'L' to a lower level nL is 
given by, 


Pn' L’ ,nL aing 

ii 

3 3. 

3 

(32) 


sing 


Pn' L' ,nLtrip 

II 

CO 

3 < 

3 

(33) 


trip 


Pn' L' ,nLtot 

Pn' L' ,nZ/sing P Pn' L' ,nLtrip 

(34) 


Pn 


L’ = T 


3 8{L')T e 


4 A. 


n ' L Li P 


(37) 


Finally the probability that either constituent particle 
will decay is 


Pn' 


n'L' 


= 2 


- V,T 


4 Ann,’ 


3 

4 A n’L' trlp ) 


(38) 


where r MiT is the decay rate of the /i or r as appropriate. 


B. Total branching ratios 

The previous subsection gave the branching ratios for 
the various decays from a particular state nL. We now 
need to consider the probability of leptonium being in 
the state nL in the first place. There are two possibili¬ 
ties for populating the state: it may recombine directly 
into the state nL , or it may recombine into a higher state 
and cascade down into nL. That is, we assume Case 
A recombination^ because leptonium is short-lived, 
there is no resonant scattering of emission lines from Ps, 
M or T, through the excitation and re-emission of sur¬ 
rounding Ps, M or T atoms. For the same reason, we 
also assume that collisional excitation of the atoms is 
negligible. 


1. Cascade probabilities 


where the leading factors of 1/4 and 3/4 are the proba¬ 
bilities that the onium will form in the singlet or triplet 
state, respectively. 


The probability of cascading from a level n'L' to a 
lower level nL via all cascade paths, C n 'L',nL , is calcu¬ 
lated using an iterative procedure!!, such that, 


J 


1, 


n = n’ A L = L' 


C , T > T = < n '- 1 

' E E P n'L',n"L" t .C„"L",n L!sias , trip , U ± n' V L ^ L ', 


(39) 


n"=n L"=L'±1 


sing,trip 

r 


where the subscripts ‘sing’ and ‘trip’ are used to specify unfortunate but necessary complication since the decay 
the calculations for the singlet state or triplet state, an 
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paths are different for each. 


Let us then define f n L as the relative population of the 
level nL , i.e., 


2. Population of the state, nL 


Meanwhile the state is depopulated by the processes 
described in section[V] Therefore, in equilibrium the pop¬ 
ulation of any state, N u l is given by[43l l47j| . 


OO 71— 1 

^ ' y ] -^pairs/ion<?rr/o- ni C n 'L',nL aiDgitrip 
n'=n L '=0 

-^^TlZ/sing^rip "^^-^singjtrip 5 


fTlL s i n 


NnL s: 


-^pairs«Aon$rr 
oo n — 1 

^ ^ v f <7 n i L r ,7lZ/ s ing,trip 


n'=n L'—O 


A„ 


(41) 


which is the fraction of leptonium in the level nL for a 
population whose formation and decay rates are in equi¬ 
librium. 


where N pa i IS is the rate at which pairs are produced, f lon 
is the fraction of pairs produced with energy less than the 
ionization energy, i.e. / 77 or f ee in sections flll Al and lHI Bl 
above, g TI is the fraction of those pairs which radiatively 
recombine (§1330 Figure [TJ, f a nL is the fraction of those 
pairs which form in level nL, 1? IIV A1 equation 1191) . and 
again the equation has two forms reflecting the different 
decay paths of the singlet and triplet states. 


3. Branching ratios 

The final step in calculating the total branching ratios 
is simply to multiply the fractional population each state, 
fnL , by the appropriate decay constant, and sum over all 
levels. Thus, 


J 


Phya — 
P Ly/3 = 
-^Baa = 


-pBa (3 — 
P ~ 


P — 

jl — 


Pee = 


decay 


A 


21,10 


A 


31,10 


— foi -I-fo-l . 

^ J -^J-sing 1 ^ J ^-*-trip J ‘ 

— H-/q 1 . 

^ J OXsing ' ^ J *- ) - L trip j ‘ 


A 


31,20 


1 


A, 


41,20 


r /3 l sinK + t /31 


^ f Using 4"~ ^ f lltrip 


A 


30,21 


3 

4 

^40,21 ( 7/4 Osing + -j; 


^ f 30 s i n g 4” ^ f 30 trip J i 


f 40 trip I ‘ 


y ] ,fn0 ain Pl( n ), 


71 — 0 
OO 


E 7 /«o trip r 3 (n), 


71 — 0 
OO 


i/nOtrip Tee (n), 


71 = 0 
OO 71—1 


71=0 L=0 


EE ( ifnL s i n g 4” a fnL 


trip I i 


Lyman a 

(42) 

Lyman /3 

(43) 

B aimer a 

(44) 

B aimer /? 

(45) 

Two photon annihilation 

(46) 

Three photon annihilation 

(47) 

Electron — positron annihilation 

(48) 

Decay of either particle 

(49) 


All of the equations 02] 09] contain either one or 

two infinite sums. All include the infinite sum in equa¬ 
tion 0T] which accounts for recombination into any level 
up n = 00 followed by cascade to the relevant level. 
Equations 0S] 0H also include the sum over all levels 
up to n = 00 which all contribute to the relevant decay 
process. In practice we cannot compute these infinite 
sums, and instead compute equations 05] 01]up to cer¬ 


tain limits on each sum. We then alter these limits, and 
fit to the resulting curve or surface, and extrapolate our 
results to n = 00 . The limits are n = 11 — 15, for the 
capture-cascade sum, and n = 6 10 for the sum over 

the decays. 

For example, Figure 0] shows the fit to the branching 
ratio for M emitting Lyman a radiation at a gas temper¬ 
ature of T = 10 4 K. We find a good fit to all curves using 
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FIG. 3: The branching ratio of Ly a emission from M at a gas 
temperature of T = 10 4 K as a function of the limiting value of 
the sum over n in equation S3 shown by the points. We find 
good fits with a function of the form given in equation 1501 
shown by the curve; in this case giving a limiting value as 
n —» oo of 0.45. 



FIG. 4: The branching ratio of decay into an electron-positron 
pair from M at a gas temperature of T = 10 4 K as a function 
of the limiting value of the sums over n in equations [41] and 
m shown by the points. We find good fits with a function 
of the form given in equation 1511 shown by the curve; in this 
case giving a limiting value as n —> oo of 0.59. 


a function of the form 

b = a + ce _3 , (50) 

where a, c, and d are all parameters to be fit. Figured] 
shows the two dimensional fit over both limits for anni¬ 
hilation into electron-positron pairs for M at T = 10 4 K. 
Again, we find a good fit with a functional form 

b = a + ce~~^ + fe ~~^, (51) 

where where a, c, d , f and g are all parameters to be fit, 
and the n\ and ri 2 are the different limiting values of n 
in the sums. 


We have checked the accuracy of the extrapolations 
in the following way. The branching ratios for the two 
photon or three photon decay of Ps, should be 0.25 and 
0.75 by definition, since there is no other way for Ps to 
decay. We find the values given in Table El which agree 
to within 2 per cent in all cases. 

We wish to highlight one further caveat to these 
branching ratios. The extrapolations for photon decays, 
and electron-positron decays are not large, and become 
less important as n increases, since they can only occur 
when L = 0, which becomes less common for higher n. 
However, the decay of either particle can occur from any 
L level, and the extrapolations are large and therefore 
uncertain. Therefore, rather than rely on the extrapola¬ 
tions for decay, we set the probability to be one minus the 
sum of the other decays. In any case, the decay of either 
particle does not have a clear observational signature and 
it is the probabilities for the other processes which are 
most important for this study. This uncertainty in our 
results could be lessened through calculations to higher 
levels, but would require significant computing resources 
which are unavailable to us at this time. 

The results for the branching ratios are given in Ta¬ 
bles 0 ED an d IVIII Since Ps is well studied and under¬ 
stood, we do not discuss it further [13. l22j|. 

M has significant probabilities of Lyman a and Balmer 
a radiation before decaying. The decay processes are 
dominated by two-photon annihilation for para-M and 
electron-positron annihilation for ortho-M, as well as sig¬ 
nificant probabilities for either ^ to decay if T < 10 5 K. 

The much shorter lifetimes of r means that the decay 
processes of T are dominated by the decay of either , 
with only a small probability of either electron-positron 
or two photon annihilation. The probabilities of emitting 
Lyman a radiation before decay are 1-2 per cent, and 
s=s 0 for other emission lines. 


VII. OBSERVABLE SIGNATURES 

Tables I VII and IVIII give the branching ratios for pos¬ 
sible radiative recombination lines and decays of M and 
T. We now discuss the observational signatures of these 
processes. The energies of radiative recombination lines 
and of 2-photon annihilation signatures and ionization 
are given in Table [T] In addition to the 2-photon anni¬ 
hilation signature there will be a continuum of radiation 
from the 3-photon decay of the singlet state, with en¬ 
ergies from zero to the energy of the 2-photon decay. 
The decay mechanisms for both r and /i are numerous 
and there is no one clear observable signature that we 
can associate with either of these, nor is there any single 
identifiable signature of annihilation into e^. 

All these signatures are observable with current in¬ 
strumentation if they are sufficiently bright. For exam¬ 
ple Fermi-LAT is sensitive to y-rays between 30 MeV 
to 300GeV and could observe the annihilation of M and 
T; INTEGRAL IBIS is sensitive to the T Lyman lines; 
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TABLE V: Total branching ratios for Ps as a function of temperature. 

Two Three 
photon photon 


T (I<) 

Lya Ly/3 Balmer a 

Balmer /3 decay 

decay 

1000 

0.46 0.08 0.24 

0.06 

0.24 

0.73 

10000 

0.29 0.07 0.09 

0.03 

0.25 

0.74 

100000 

0.11 0.03 0.03 

0.01 

0.25 

0.75 

1000000 

0.04 0.01 0.02 

0.01 

0.25 

0.75 

10000000 

0.03 0.01 0.02 

0.01 

0.25 

0.75 


TABLE VI: Total branching ratios for M as a function of temperature. 

Two Three 

photon photon Electron-positron Decay of 

T (K) Lya Ly/3 Balmer a Balmer j3 decay decay decay either particle 


1000 

0.37 0.04 0.27 

0.04 

0.16 

0.00 

0.48 

0.36 

10000 

0.45 0.05 0.31 

0.05 

0.20 

0.00 

0.59 

0.21 

100000 

0.46 0.07 0.26 

0.06 

0.23 

0.00 

0.70 

0.07 

1000000 

0.32 0.07 0.12 

0.04 

0.25 

0.00 

0.74 

0.01 

10000000 

0.13 0.04 0.02 

0.01 

0.25 

0.00 

0.74 

0.00 


XMM-Newton is sensitive to X-rays between 0.1 and 15 
keV, and can observe the recombination lines of M, and 
the T Balmer lines. 


VIII. EXPECTED SIGNATURES FROM 
ASTROPHYSICAL SOURCES 

We are finally in a position to estimate the expected 
observable signatures of M and T from astrophysical 
sources. Consider a source of particle anti-particle pairs 
being produced at a rate r per second. The fraction of 
these pairs which will produce M or T is given by 

,/oniimi — C./ion Qrr (52) 

where /i on is the fraction of pairs produced with energy 
less than the ionization energy, i.e. / 77 or f ee in sec¬ 
tions IIII Al and IIII Bl above, and g TI is the fraction of 
those pairs which radiatively recombine (§ Ml Figure [lj . 
Therefore the luminosity of any particular observable sig¬ 
nature in ph s -1 is 

L — P /onium b (53) 

where b is the branching ratio of the particular signature 
(Tables IVTl and [Ylll) . and p is a factor which specifies the 
number of photons produced by the process, i.e. p = 1 for 
recombination lines, p = 2 for two-photon annihilation 
and p = 3 for three-photon annihilation. These lumi¬ 
nosities can be easily converted to units of erg s -1 using 
the energies of the emitted photons given in Table HI 
We will now estimate the expected fluxes from partic¬ 
ular sources. These should be compared to the limiting 
sensitivities of current instruments. The M recombina¬ 
tion lines are all in the soft X-ray band (0.5 - 2 keV), 
for which a 100 ksec observation has a ss 4cr limiting 
sensitivity (48| of /x ~ 3.1 x 10 “ 16 erg cm -2 s -1 . The T 


recombination Balmer lines are in the hard X-ray band 
(2 - 10 keV), for which a 100 ksec observation has a 
ss 4 a limiting sensitivity of /x ~ 1.4 x 10 “ 15 erg cm “ 2 
s -1 . The T Lyman lines are at the low energy limit of the 
SPI spectrometer on board INTEGRAL. The limiting do- 
sensitivity for an exposure time of 10 5 s is « 3 x 10 “ 13 
erg cm -2 s -1 , using data from the INTEGRAL web¬ 
pages. The M and T annihilation lines are detectable 
by FERMI-LAT, for which the limiting 4 a sensitivity for 
an exposure time of 10 5 s is ~ 4 x 10 “ 13 erg cm -2 s _1 
for M and « 1 x 10 “ 13 erg cm -2 s _1 for T, using data 
from the FERMI webpages. 


A. Jets 

The powerful relativistic jets of active galactic nu¬ 
clei (AGN) and microquasars are a probable source of 
pair production. These jets have strong radio emission 
with a power-law spectrum which is interpreted as being 
due to synchrotron emission from relativistic e~ gyrat¬ 
ing around the magnetic field lines of the AGN. The jets 
must be electrically neutral, otherwise they would cause 
a potential difference to build-up, which would oppose 
and eventually stop the jet. 

It is unknown whether the positive component of jets 
consists of positrons, protons, or a mixture of both. Ob¬ 
servational studies have had to rely on indirect methods 
of searching for the presence of e + , such as estimates 
of the bulk kinetic energy contained in jets, which have 
been used to argue for both e“-p plasma |49| and e“-e + 
plasma (50l - l52| . 

Theoretically there are good reasons to expect that 
jets contain some fraction of c + . A pair plasma has the 
advantage of explaining 7 -ray jets[53j and the very high 
Lorentz factors (F > 5) required to account for superlu¬ 
minal bulk velocities of jetsjU- 









12 


TABLE VII: Total branching ratios for T as a function of temperature. 

Two Three 

photon photon Electron-positron Decay of 

T (K) Lya Ly/3 Balmer a Balmer /3 decay decay decay either particle 


1000 

0.01 0.00 0.00 

0.00 

0.03 

0.00 

0.06 

0.91 

10000 

0.01 0.00 0.00 

0.00 

0.03 

0.00 

0.07 

0.90 

100000 

0.01 0.00 0.00 

0.00 

0.04 

0.00 

0.09 

0.87 

1000000 

0.01 0.00 0.00 

0.00 

0.06 

0.00 

0.12 

0.83 

10000000 

0.02 0.00 0.00 

0.00 

0.08 

0.00 

0.17 

0.75 


It is possible that some Ps may form in the jet, or 
similarly, that some e^ pairs may collide with the re¬ 
quired energy to form M and T as outlined in § IIII Bl 
However, in the observers frame, the velocity of the oni- 
ums thus formed will be highly relativistic, causing rela¬ 
tivistic Doppler broadening and beaming of any emitted 
radiation [55], [56j. Such a broadened signal may still be 
possible to observe[55|, but would be difficult to inter¬ 
pret, and hardly constitutes an unambiguous test of the 
presence of exotic oniurn atoms such as M or T. There¬ 
fore we do not consider onium formation in the jet itself 
as a likely candidate for detection. 

However, if a jet collides with stationary object such 
as a gas cloud, or a star, then the leptonium formed in 
the collision will give rise to emission that is in principle 
observable. Such collisions may occur in the radio jets of 
AGNptfl, or if the jets of microquasars are mis-aligned 
and hit the secondary companion|58|. 

Let us consider some illustrative examples. 
Previously^, we have calculated the expected positron 
contents of jets, by applying the arguments of Marscher 
et al. [56j] and Marscher fSQj , which were developed in 
order to search for 511 keV Ps annihilation radiation 
from 3C 120, to the empirical measurements of Ghisellini 
et al. [gcj. We found a maximum of ps 10 49 e + s _1 
are produced in the jets of blazars, whereas quasars 
produce « 10 46 - 10 48 e + s _1 in their jets. The jets of 
microquasars are expected to have a positron flux of 
ps 10 41 e + s _1 f58|. We assume that the object being hit 
by the jet is dense enough to stop all these positrons. 
Let us assume a spectral index of a = 1.5 and a temper¬ 
ature of T = 10 6 K for the positrons in the jet. Thus 
the intrinsic luminosities of the various observational 
signatures can be calculated, and are given in Table IyTh! 
for a blazar. These are then used to calculate the flux of 
the signal as a function of redsliift, as shown in Figure [5] 
Similar results are shown in Table HXl and Figure [Gj for a 
mis-aligned microquasar with a spectral index of a = 1.5 
and a temperature of T = 10 5 K for the positrons in the 
jet. 

For our illustrative examples, no lines from blazars are 
above the relevant current 4cr detection limits for a 100 ks 
exposure for either M or T. However, a microquasar jet 
which produces 10 41 e + s _1 , with spectral index= —1.5 
and T = 10 5 K, colliding with a gas cloud, would be 
detectable in principle to the distances given in Table [X] 
The recombination lines would only be detectable at very 



(a)M 



z 

(b)T 


FIG. 5: The flux as a function of redsliift for a blazar jet 
which produces 10 49 e + s _1 , with spectral index a = 1.5 and 
T = 10 6 K, colliding with a gas cloud. The lines show two- 
photon annihilation (red), Ly a (black) Ly /3 (dashed, black), 
Ba a (blue) and Ba /3 (dashed, blue). 
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TABLE VIII: The intrinsic luminosities for a blazar jet which 
produces 10 49 e + s _1 , with spectral index a = 1.5 and T = 
10 6 K. 

Luminosity (ph s _1 ) 



M 



T 



Ly a 

2.7 

X 

hP" 

1.2 

X 

Tip 

Ly P 

6.0 

X 

10 41 

1.3 

X 

io 40 

Ba a 

9.8 

X 

10 41 

7.7 

X 

10 39 

Ba p 

3.3 

X 

10 41 

1.6 

X 

10 39 

Two photon 4.1 

X 

10 42 

9.6 

X 

10 41 


TABLE IX: The intrinsic luminosities for a microquasar jet 
which produces 10 41 e + s _1 , with spectral index a = 1.5 and 
T = 10 5 K. 


Luminosity (ph s 4 ) 



M 



T 



Ly a 

3.9 

X 

TP 

0° 

X 

TP 

Ly P 

5.8 

X 

10 33 

9.7 

X 

10 31 

Ba a 

2.3 

X 

10 34 

5.9 

X 

10 31 

Ba P 

4.8 

X 

10 33 

1.2 

X 

10 32 

Two photon 4.0 

X 

10 34 

7.2 

X 

10 33 


close distances, but the two photon annihilation should 
be detectable at 4er in a 100 ksec exposure to « 0.4 kpc 
and « 1 kpc for M and T respectively. 


B. Accretion discs 

We now consider pairs produced in the accretion disc 
itself through photon-photon annihilation. The number 
density of pairs thus produced in the optically thick re¬ 
gion of the disc is 


7V 77 


1 

or R 


(54) 


where <tt is the Thompson cross-section for /i or r and 
R is the radius of the optically thick disc(filj|. If we as¬ 
sume R « 2 GM/c 2 , i.e. the disc has approximately the 
Schwarzschild radius, and the thickness of the disc is « f?, 
then taking masses of Magn = 10 6 M 0 for the mass of a 
black hole in an AGN, and qso = 10 M 0 for the mass 
of a black hole in a microquasar, then the pair yields are 
as given in Table |XH 



kpc 

(a)M 



kpc 

(b)T 

FIG. 6: The flux as a function of redshift for a microquasar 
jet which produces 10 41 e + s _1 , with spectral index= —1.5 
and T = 10 5 K, colliding with a gas cloud. The lines show 
two-photon annihilation (red), Ly a (black) Ly p (dashed, 
black), Ba a (blue) and Ba p (dashed, blue). 


TABLE X: The distance to which a microquasar jet which 
produces 10 41 e + s _1 , with spectral index= —1.5 and T = 
10 5 K, colliding with a gas cloud, would be detectable at 4 a 
in 100 ksec. 

Distance (kpc) 

M T 

Ly q 0.04 9 x 10~ 4 

Ly /3 0.01 3 X 10" 4 

Ba a 0.01 0.001 

Ba j3 0.007 7 x 10 -4 

Two photon 0.40 1.13 


Now, assuming a spectral index of a = 1.5, then we 
may estimate the branching ratios as in equation 1531 
whereupon we find the luminosities in ph s , as listed in 
Table IXTT1 The flux of the AGN as a function of redshift, 
and the flux of the microquasar as a function of distance, 
are shown in Figures Q] and Figures [8] 

For these illustrative examples, and a 4cr detection 
limit for a 100 ks exposure, the two-photon annihilation 
line of M would be observable in the accretion disc of 
an AGN at z < 0.025, while the two-photon annihilation 
line of T would be observable out to z < 0.9. Neither the 
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TABLE XI: The number of pairs produced through photon- 
photon annihilation in an accretion disc. 

AGN Microquasar 
/t 1.7 x 10"- 1.7 x ID 4 - 

t 5.0 x 10 54 5.0 x 10 44 


TABLE XII: The intrinsic luminosities for pair production 
via photon annihilation in accretion discs surrounding a mi¬ 
croquasar and an AGN, assuming solar masses of 10 and 10 6 
respectively, temperatures of 10 5 K and 10 6 K, and a 
spectral index of a = 1.5. 

Microquasar AGN 

Luminosity (ph s _1 ) Luminosity (ph s _1 ) 



M 



T 



M 



T 



Ly a 

2.9 

X 

10 3t> 

1.8 

X 

10 3B 

2.0 

X 

IF 

2.4 

X 

TcF 

Ly d 

4.2 

X 

10 34 

2.0 

X 

10 35 

4.4 

X 

10 44 

2.7 

X 

10 45 

Ba a 

1.7 

X 

10 35 

1.2 

X 

10 35 

7.2 

X 

10 44 

1.6 

X 

10 45 

Ba /3 

3.6 

X 

10 34 

2.5 

X 

10 34 

2.4 

X 

10 44 

3.2 

X 

10 44 

Two photon 2.9 

X 

10 35 

1.5 

X 

10 37 

3.0 

X 

10 45 

2.0 

X 

10 47 


M nor the T recombination lines would be detectable. 

Meanwhile, for microquasars the M and T lines would 
be detectable at 4cr in a 100 ksec exposure to the dis¬ 
tances given in Table IXIIII The recombination lines are 
detectable at close distances while the two photon anni¬ 
hilation should be detectable to 1.1 or 51 kpc for M and 
T respectively. 


IX. CONCLUSIONS 

True muonium and true tauonium are the most com¬ 
pact pure QED systems, but have never been observed. 
Unlike Ps, for which there are extensive observations in 
our own Galaxy [T3|. both the formation and decay of M 
and T are affected by the intrinsic instability of the /r 
and r leptons. We have investigated the likelihood of 
their formation in astrophysical environments, and the 
prospects for their observation. 

The probability of formation is small, ~ 10 - ' from 
photon-photon annihilation or electron-positron annihi¬ 
lation. The probability is small for two reasons: (i) the 
lifetimes of the n and r are intrinsically short (S IVBI) . and 
thus M and T can only form from the products of pair 


TABLE XIII: The distance to which an accretion disc around 
a microquasar with a 10 Mg black hole, and spectral index 
a = 1.5 and T = 10 s K, would be detectable at 4 a in 100 
ksec. 

Distance (kpc) 

M T 


Ly a 

0.11 

0.04 

Ly /3 

0.05 

0.01 

Ba a 

0.04 

0.06 

Ba (3 

0.02 

0.03 

Two photon 

1.1 

51 



z 

(a)M 



z 

(b)T 


FIG. 7: The flux as a function of redshift for an accretion 
disc around an AGN with a 10 6 Mg black hole, and spectral 
index a = 1.5 and T = 10 6 K. The lines show two-photon 
annihilation (red), Ly a (black) Ly /3 (dashed, black), Ba a 
(blue) and Ba /? (dashed, blue). 


production processes, such that the pairs immediately 
recombine, (ii) even then, only those pairs with a total 
kinetic energy less than the ionisation energy can form 
leptonium, and since pair production usually takes place 
in high energy process, these pairs constitute a small frac¬ 
tion of the total. Nevertheless, high energy astrophysical 
environments are capable of producing copious numbers 
of /i and r pairs, and the cross-section for radiative re¬ 
combination dominates that of direct annihilation (§EI- 
Thus even the small fraction of pairs with energy low 
enough to form M or T can lead to a significant flux. 

The decay of M and T are hastened by the short life- 
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kpc 

(a)M 



(b)T 

FIG. 8 : The flux as a function of redshift for an accretion disc 
around a microquasar with a 10 M© black hole, and spectral 
index a = 1.5 and T = 10 5 K. The lines show two-photon 
annihilation (red), Ly a (black) Ly /3 (dashed, black), Ba a 
(blue) and Ba /3 (dashed, blue). 


times of the /i and r leptons, which may be one rea¬ 
son why the possibility of astrophysical observations has 
not received much attention. Here, we have carefully 
calculated the probabilities of the various observational 
signatures. We have calculated the cross-sections for re¬ 
combination onto the nLth level (§ II V All , and the re¬ 


sulting branching ratios for Lyman a, Lyman /3, Balmer 
a, Balmer /3, two photon annihilation, three photon an¬ 
nihilation, annihilation into e 1 * 1 pairs, and the decay of 
either /( or r. Although the decay of either r does domi¬ 
nate for T, there is still a small probability for observing 
recombination lines or two-photon annihilation. For M 
the situation is more hopeful, with significant branching 
ratios for Lyman a, Balmer a and two-photon annihila¬ 
tion. 

In section IVIIII we made estimates of the fluxes of M 
and T recombination and annihilation signatures for the 
cases of blazar jet-cloud interactions, jet-star interactions 
in mis-aligned microquasars, and within the accretion 
discs of AGN and microquasars. These were compared to 
the current detection limits of X-ray and 7 -ray observato¬ 
ries. The expected signatures from AGN jet-cloud inter¬ 
actions were all below current detections limits. However, 
M and T formation within microquasar jet-star interac¬ 
tions, or within the accretion discs of both AGN and 
microquasars were estimated to yield signatures brighter 
than the detection limits, with those from microquasars 
offering the brightest estimates due to their proximity. 
Actual observations would be further complicated by the 
intrinsic backgrounds and also emission from the object 
in question, which may have other significant line emis¬ 
sion (e.g. see Marshall et al. [621] for a spectrum of SS 
433). 

These examples are only illustrative. Other sources 
may also be significant. For example, 511 keV electron- 
positron annihilation radiation was recently discovered 
in terrestrial 7 -ray flashes due to lightning strikes Q. 
These same flashes could lead to M or T formation, 
as well as the possibility of observing Ps recombination 
lines. As a different example we note in passing that 
the recent detection of an unidentified line at ss 3.5 keV 
in galaxy clusters [H, [3 cannot be explained by the T 
Balmer a line at 3.3 keV, which is ruled out by the con¬ 
straints on the energy. 

In summary, the astrophysics of pair-production in any 
high energy source could lead to possible M and T for¬ 
mation. This paper provides the tools to estimate the 
fluxes of the M and T recombination and annihilation 
signatures once the rate of pair production is known. 
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Appendix: Energy threshold in the laboratory frame for muon pair production from electron-positron 

annihilation 

To relate the lab frame energies of the colliding electrons, E\ and E 2 , to the zero momentum frame Lorentz factors, 
7, of the produced muons, we note that the absolute square of the relativistic four momentum p 2 = p^p^ is invariant. 
In the lab frame, let the incoming electron be traveling along the x-axis in the positive direction, with speed u, and 
the incoming positron be at an angle 0 to the x-axis, in the x-y plane, with speed w. The relative speed of the particles 
in the lab frame is then, v = u — w cos 0. Therefore in the lab frame the total relativistic four-momentum of both 
particles is given by, 


Plab 


Ei + E2 


, pi +p 2 cos 0, p 2 sin 0,0 I . 


(A.l) 


In the zero-momentum frame, after the collision, the produced muons have equal and opposite momentum, thus, 

£3 + £4 


Pz.ra. — 


-, 0 , 0,0 . 


(A.2) 


Now since p 2 is invariant we have, 


9 2 

P\ab Pz.m. 


Ei + E2 


— {pi + P 2 cos 0) 2 — p\ sin 2 0 = 


£3 + £4 


£1 £ 2 + c 4 (?n 2 — 2M 2 y 2 ) = \J~Ef- cos0. (A. 3) 


Therefore the relative velocity is, 


sjc 2 E 2 — m 2 c 6 \Jc 2 E\ — m 2 c 6 cos 0 


£1 


£2 


(A.4) 
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